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The origin of pulsar kicks is reviewed in the framework of the spin-flip conversion of neutrinos 
propagating in the gravitational field of a magnetized protoneutron star. We flnd that for a mass 
in rotation with angular velocity uj, the spin connections entering in the Dirac equation give rise 
to the coupling term u: ■ p, being p the neutrino momentum. Such a coupling can be responsible 
of pulsar kicks owing to the neutrino emission asymmetry generated by the relative orientation of 
p with respect to uj. For our estimations, the large non standard neutrino magnetic momentum, 
/i^ < 10~^^/i_B, is considered. 
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I. INTRODUCTION 



' The origin of the pulsar velocity represents till now an unsolved and still discussed issue of the modern astrophysics. 
Observations show that pulsars have a very high proper motion (velocity) with respect to the surrounding stars, with 
a three-dimensional galactic speed of 450±90Km/sec up to values greater than fOOOKm/sec [1-6]. This suggests that 
' nascent pulsars undergo to some kind of impulse (kick). After the supernova collapse of a massive star, neutrinos 
carry away almost all (99%) the gravitational binding energy (3 x lO^^erg). The momentum taken by them is about 
lO'^^gr cm/sec. An anisotropy of ~ 1% of the momenta distribution of the outgoing neutrinos would suffice to account 
fS| . for the neutron star recoil of 300Km/sec. Even though many mechanisms have been proposed, the origin of the pulsar 
T— I ' kick remains an open issue. 

An interesting and elegant mechanism to generate the pulsar velocity, which relies on the neutrino oscillation physics 
in presence of an intense magnetic field, has been proposed by Kusenko and Segre (KS) [7]. The basic idea is the 
I following. Electron neutrinos Vf. are emitted from a surface placed at a distance from the center of a pulsar greater 
' than the surfaces corresponding to muon neutrinos Vfj, and tau neutrinos Vr . Such surfaces are called neutrinospheres. 
Oh^ Under suitable conditions, a resonant oscillation j/g — » Vf^^r may occur between the Ve and v^^^r neutrinospheres. 
Q . The neutrinos are trapped by the medium (due to neutral and charged interactions) but neutrinos f^,r generate 
?H ' via oscillations can escape from the protoneutron star being outside of their neutrinosphere. Thus, the "surface 
c/j . of the resonance" acts as an "effective muon/tau neutrinosphere". In the presence of a magnetic field (or some 
' other nonisotropic effect), the surface of resonance can be distorted and the energy flux turns out to be generated 
^ . anisotropically. In [7] , the anisotropy of the neutrino emission is driven by the polarization of the medium due to the 
JJ"! magnetic field B. The usual MSW resonance conditions turn out to be modified by the term [8-12] 



where p = p/p, p is the neutrino momentum, e is the electric charge, Gj? is the Fermi constant, and n,, is the electron 
density. Nevertheless, the neutrino masses required for the KS mechanism seem to be inconsistent with the present 
limits on the masses of standard electroweak neutrinos. These limits do not apply to sterile neutrinos (they may have 
only a small mixing angle with the ordinary neutrinos) [13,14]. Papers dealing with the origin of pulsar kicks can be 
found in [15-39]. 

The aim of this paper is to suggest a further mechanism to generate pulsar kicks. It is based on the spin flavor 
conversion of neutrinos propagating in a gravitational field generated by a rotating source. Even though the grav- 
itational field per se cannot induce neutrino oscillations, unless a violation of the equivalence principle is invoked 
[40,41], it affects the resonance conditions (hence the probability that left-handed neutrinos convert into right handed 
neutrinos, the latter being sterile may escape from the protoneutron star). Such a modification is induced by spin 
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connections entering in the Dirac equation in curved spacetimcs. They give rise to the couphng term ^ u; • p, where 
U3 is the angular velocity of the gravitational source. The relative orientation of neutrino momenta with respect to 
the angular velocity determines an asymmetry of the neutrino emission, hence it may generate the pulsar kicks. 

The paper is organized as follows. In Sect. 2 we shortly review the Dirac equation in curved space-times. Sect. 3 
is devoted to the computation of the fractional asymmetry. The resonance and adiabatic conditions, as well as the 
spin flip probability are also discussed. Conclusions are drawn in Sect. 4. 



II. DIRAC EQUATION IN THE LENSE-THIRRING GEOMETRY 

The phase of neutrinos propagating in a curved background is generalized as [42-52] 

3 

where / is the flavor index and j the mass index. U fj are the matrix elements transforming flavor and mass bases. 
P ■ Pnuii = PfiPnuii > where is the four-momentum operator generating space-time translation of the eigenstates 
and p^^;; = dxf^/dX is the tangent vector to the neutrino world-line x^, parameterized by A. The covariant Dirac 
equation in curved space-time is (in natural units) [53] [i"/f^{x)Df^ — m]'ip = 0, where the matrices 7^(2;) are related 
to the usual Dirac matrices 7" by means of the vierbein flelds e^(a;), i.e. ^^{x) = e^(a;)7°. The Greek (Latin 
with hat) indices refer to curved (flat) space-time. is defined as D^j, = + r^(a;), where T^{x) are the spin 

connections r^(a;) = \['y°' (semicolon represents the covariant derivative). After some manipulations, the 
spin connections can be cast in the form [42] 

7''r^ = 7"e^ {iAc, [-{-9)-"' l"] } , (2) 

where 

A'g = \^9ele''Hei^,^ - , (3) 

g = det{gij,^). For geometries with a spherical symmetry, such as the Schwarzschild space-time, Aq vanishes. For 
rotating gravitational sources, instead, Aq acquires non vanishing components. To prove this, let us consider the 
Lense-Thirring geometry, i.e. the geometry of a rotating source. The line element (in the weak field approximation) 

is given by 

ds^ = (1 - (p){dtf -{1 + (/>)(dx)2 - 2h • dxdt , (4) 

where 

2GM AGMR'^ 
x = {x,y,z), (f>=——, h= a; Ax, 

M is the gravitational mass, and R its radius. The vierbein fields associated to the metric (4) are 

eoo = (l + '^)' e.Qi = -h\ €■■= -{I - <j))5ji , 



Prom Eq. (3) one finds 



where 



, 3(a) • x) x 

U} =UJ TT-^ 



The operator 7^r^ ~ j'^AcfiJ^ acts differently on left- and right-handed neutrino states [42]. By writing 7^ = 
Vr — Vl, where Vl.r = (lT7^)/2, one immediately sees that left- and right-handed neutrinos acquire a gravitational 
potential which is opposite for the two helicities. The dispersion relation of neutrinos turns out to be modified by 
the term [42] p^AcfiJ^ = — p • A.gJ^- In the case of neutrino oscillations, it is convenient to add, without physical 
consequences, a term proportional to the identity matrix, so that 7^ can be replaced by the left-handed projection 
operator Vl- Spin-gravity coupling terms can be then pushed in the left-handed sector of the neutrinos effective 
Hamiltonian (see Eqs. (5)- (8)). 
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III. THE ASYMMETRIC EMISSION OF NEUTRINOS FROM PROTONEUTRON STARS 



Neutrinos inside their neutrinospheres are trapped by the weak interactions with the background matter, 
therefore acquire the potential energy 



They 



*gr cm 



■y/(r, i)eV, 



where / = e,/i,T, p is the matter density, yg =Y^ — 1/3 and y^^T = — 1. In these expressions, Yg is the electron 
fraction. In the present paper we shall consider the case in which matter induced effective potential is \V^J <C 1 (as 
shown in [54], Vi^^ may cross the zero at r ~ 12km). This occurs in the regions where the electron fraction Yg assumes 
the value 1/3 [ye < 1) [55-60,39,54]. 

Protoneutron stars posses large magnetic fields which may vary in a range of several order of magnitude. In such 
astrophysical systems, the interaction of neutrinos with (uniform) magnetic fields is described by the Lagrangian 
[61-63] 

ab 1 



where jj, is the magnetic momentum of the neutrino, -F-g is the electro-magnetic field tensor, and a' 



db 



lb", 7^ 



Spin flavor conversions of neutrinos propagating in magnetic fields are generated by the transition magnetic moments, 
which are non-diagonal terms in the effective Hamiltonian describing the neutrino evolution. 

Taking into account the gravitational and magnetic interactions, the equation of evolution describing the conversion 
between two neutrino flavors / = e and /' = /i, r reads (we refer to Dirac neutrinos, but the formalism also applies 
to Majorana neutrinos) [64] 



dX 



( "f^ \ 

Uf'L 

\ Vf'R I 



n 




(5) 



where, in the chiral base, the matrix Ti is the effective Hamiltonian 



n 



Hl n*jj, 

Tiff Hr 



(6) 



K. 



-no 

6S2 



SC2 6S2 

K^, + + 6c2 



(7) 



— 6C2 6S2 
SS2 6C2 



IJ-ff fj-ff 



Pi.A'^ir) _ AGMR^ 
El ~ 5r^Ei 

13 / R y/lOkmVcycos^ 



p • U) 

2 



(8) 



(9) 



Here 5 = (Am^ = m| — mf), C2 = cos 2^, 82 = sin 20, 9 is the vacuum mixing angle, Ei is the energy measured 
in the local ^ame, B_l = B sin a is the component of the magnetic field orthogonal to the neutrino momentum, and 
(3 is the angle between the neutrino momentum and the angular velocity. 

r^G is diagonal in spin space, so that it cannot induce spin-flips. Its relevance comes from the fact that it modifies 
the resonance conditions (resonances are governed by the 2x2 submatrix in (6) for each pairs of states). For the 
transition ffL^ ^S'R ^"^^ ™ ^^^^ §®*^ 



Ke+fiG(r)-25c2 = 0, 



(10) 



where r is the resonance point. 
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Qg in (10) distorts the surface of resonance owing to the relative orientation of the neutrino momentum with respect 
to the angular velocity. As a consequence, the outgoing energy flux results modified. To estimate the anisotropy of the 
outgoing neutrinos, one needs to evaluate the energy flux Fg emitted by nascent stars. According to the Barkovich, 
D'Olivo, Montemayor and Zanella paper [65], the neutrino momentum asymmetry is given by 

^ 1 /o" F ■ u,da _lg 

The factor 1/3 comes from the fact that only one neutrinos species is responsible for the anisotropy (thus it carries 
out only 1/3 of the total energy). Fg is the outgoing neutrino flux through the element area da of the emission 
surface, w' = ui' /\u}'\ is the direction of the vector a;', whereas fi is the unity vector orthogonal to da (notice that the 
area element in curved spacetime is da; in the weak field approximation one has da ^ da) . g is the radial 
deformation of the effective surface of resonance. It shifts the resonance point f to r(0) = f+gcoscj), with g <^f and 
cos (p = u}' - p. To determine f , one uses the resonance condition 

2^C2-Ke(r) = 0. (12) 

The deformation g is evaluated by expanding (10) about to r{<p), and using the shifts p — > p + (5p, -^Vv^+ Sy^^ 
[65], where 

6p = ^ ^ p g cos (f) = hp^p g cos (j) , (13) 

Sv^^ = ^^^^"^ e cos = hy^^ gcoscj), (14) 



Eqs. (12), (13), and (14), allow to write Eq. (11) as 

Ap AGMB?uj' 



p 5r3 W^j{hp^ + hyl ) 



(15) 



To compute + hy^ , a model for the protoneutron star has to be specified. For our purpose, we adopt the 
polytropic model. The inner core of a protoneutron star is consistently described by a polytropic gas of relativistic 
nucleon with adiabatic index F = 4/3 [66]. As in [65], we assume that such a model also holds for the rest of the star. 
The pressure P and matter density p are related by P = Kp^ [66,65], where K = Tc/rrinpV^ — 5.6 x 10~^MeV~^/^. 
Tc = 40MeV and pc — IQ^^gi/cTc? are the temperature and the matter density of the core, respectively. The matter 
density profile p{r) can be chosen as [65] 

p^-^{x)= pl-\ax^ +bx + c], (16) 

where x = Tc/r, 

a = (1 — /x)Ar , b = (2/i — l)Ar , c = 1 — p.Xr , 

and Ar = GMdT — l)/rcP^^^KT ~ 0.87. rc = 10km is the core radius and Mc c^: Mq is the mass of the core {Mq is 
the solar mass). The parameter /x is determined by setting p{Rs) = 0, 



Xr{Rs-ra) Rs 



R 



(17) 



The temperature profile T(r) is related to the matter density p{r) by the relation [66,65] 

dT^ 9kL, 

-1— = rP^ (18) 

dr 7rr^ 

where Lc ~ 9.5 x lO^^erg/sec is the core luminosity, and k ~ 6.2 x 10~^^eV~^. Eq. (18) can be integrate via (16) 



T{r) = Te V^Xcixix) - X(l)] + 1 , (19) 

where Ac = 9k Lc Pc/'^.TrT^rc ~ 1.95, and 

3 

X{x)=c^x+-bc'^x'^+c{ac + b'^)x^+ (20) 
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Assuming the thermal equilibrium between neutrinos and the medium, so that p ^ T [65], and being V^^ ~ p, one 
can rewrite the inverse characteristic lengths and hy^ as = h^^ and hy^ = h~^. Eqs. (16) and (19) imply 
(at the resonance) 

d In T p{f) f Tc X 



As before pointed out, V^^ ^ 1 {ue <C 1) at r ~ 12km. For resonances occurring at f ~ 12km, where p{f) ~ 
10"gr/cm3, Eqs. (16) and''(17) give ~ 12.5km. Prom Eqs. (19), (21), (22) and (15), one infers 

' " 6, (23) 



lO^Ve lO^Hz 

where Ap/p ~ 0.01 and u' = 2u> (i.e. w || x) have been used. Typical angular velocities of pulsars w ~ (1 — 10^)Hz 
lead to ^ 10"^** — 10~^°, as follows from Eq. (23). For neutrinos with momentum p ^ lOMeV, the resonance 
condition (12) can be recast in the form Am^ cos 26* ~ 7.6 lO-^i/eeV^, i.e. 10"^° < Am^ cos26'/eV^ < 10-*^. The latter 
is consistent with the best fit of solar neutrinos Am|„„ (10~^ -j- 10~^)eV^ and 0.8 < sin^ 26sun ^ 1 [67]. 
Some comments are in order: 

• In addition to the level crossing (10), it must also occur that the crossing is adiabatic, i.e. the adiabatic 
parameter 7, which quantifies the magnitude of the off-diagonal elements with respect to the diagonal ones of 
(6) in the instantaneous eigenstates, must be greater than one when valuated at the resonance f 12km. To 
show this, let us define 1) the precession length 

L = 



^{2iiff,Bj_Y + (K^ + - 25c2Y ' 
that at the resonance, it reduces to 

^lff'B± p,ff> B± 

and 2) the width of the resonant spin flavor precession 

= 2AA , 

being A = n/Lres^ and 



P{r)) VIMY 

A is derived assuming 2/j(r, t) = (FJ <C p' j p) [54]. The spin flavor conversion is adiabatic provided A,. Lres-, 
or equivalently 

where 

Since the magnetic fields inside to protoneutron stars arc greater than B > lO^^G, whereas pffi < 10~^^pB, as 
provided by astrophysical and cosmological constraints [68,69], one gets 7 1, i.e. the adiabatic condition is 
fulfilled (such a result follows by using Eq. (16)). 
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• The conversion probability Pj^^^. is given by [61] 

P.,^^.^,^ = 1/2 - (1/2 - P) cos2^iCOs2^/ , 
where P = e"''"^/^ is the Landau-Zencr probability, and 9 is defined as 

tan2^(r) = 2^iff,B^/{nG{r) + - 2602) . 

9i = 9{ri) is the initial mixing angle of neutrinos produced at r^, and 9f = 9{rf) is the mixing angle of neutrinos 
where the spin flip probability is evaluated. Since 7^1, the Landau-Zener probability P vanishes. 

• The weak field approximation is fulfilled since {AGMR^ /5f^)ujRs < 10"^ as < IQ-^Hz and Rs > 12 km. This 
means that the geometry of the rotating sources is correctly described by (4) up to angular velocities u) lO^Hz. 
Since the typical angular velocity of the protoneutron star varies in the range (10^ — 10^)Hz, the weak field 
approximation here used can be applied. 

• The rotational effects are relevant during the time scale to ^ 10 sec {to is the time scale for the emission of the 
energy ~ 0.5 x 10^^ erg by each neutrinos degree of freedom with p ^ lOMeV) [69] . 



IV. CONCLUSION 



In this paper it has been suggested a mechanism for the generation of pulsar kicks which accounts for the spin- 
gravity coupling of neutrinos propagating in a gravitational field of a rotating nascent star. Owing to the relative 
orientation of neutrino momenta with respect to the direction of the angular velocity, the energy flux turns out to 
be generated anisotropically. Results imply a correlation between the motion of pulsars and their angular velocities. 
Such a connection seems to be corroborated by recent statistical analysis and observations discussed in [70,71]. 

Spin-gravity coupling is strictly related to the so called gravito-magnctic effect, an effect predicted by General 
Relativity [72], as well as by many metric theories [73-75]. Its origin is due to the mass-energy currents (moving or 
rotating matter contributes to the gravitational field, in analogy to the magnetic field of moving charges or magnetic 
dipole). Experiments involving the technology of laser ranged satellites [76,77] are at the moment the favorite candidate 
to test gravito-magnetic effects. In connections with the mechanism proposed in this paper, a direct evidence of the 
gravito-magnetic effect might be provided by pulsar kicks. Future investigations on the velocity distribution of pulsars 
will certainly allow to clarify this still open issue. 

Results of this paper (as well as the papers [19,22,36], in which pulsar kicks are discussed in relation to gravitational 
waves) have been obtained in semiclassical approximation, i.e. the gravitational field is described by the classical field 
equations of General Relativity. It will be of interest to investigate within the framework of quantum gravity theories. 
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Mosquera Cuesta and P.B. Pal for their comments. Many thanks to M. Barkovich, J.C. D'Olivo and R. Montemayor 
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